Efficient Selfconsistent Calculations of Multiband Superconductivity in UPd2Al3 
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An efficient physically motivated computational approach to multiband superconductivity is in- 
troduced and applied to the study of the gap symmetry in a heavy- fermion, UPd2Al3. Using realistic 
pairing potentials and accurate energy bands that are computed within density functional theory, 
self-consistent calculations demonstrate that the only accessible superconducting gap with nodes 
exhibits d-wave symmetry in the A\g representation of the Dgh point group. Our results suggest 
that in a superconductor with gap nodes the prevailing gap symmetry is dictated by the constraint 
that nodes must be as far as possible from high-density areas. 

PACS numbers: PACS numbers: 74.25. Jb, 74.20.-z 
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Superconductivity (SC) in heavy-fcrmion (HF) mate- 
rials exhibits a fascinating complexity of phenomena [Q . 
The nature of the SC state is mostly anticipated to be 
unconventional, i.e., an additional symmetry is broken 
in the SC state, and the order parameter is therefore 
not the conventional spin singlet s-wave type ^ . The 
identification of the order parameter symmetry is an es- 
sential issue in current studies of unconventional SC. A 
common approach to address the symmetry of the or- 
der parameter is through investigations of the SC gap, 
which is expected to have the same symmetry as the or- 
der parameter. While progress has been made for high- 
Tc superconductors (HTSC), the identification of the gap 
symmetry remains an open question for most HF mate- 
rials. The eventual identification of the gap symmetry 
in a HF compound is, however, not sufficient to identify 
the pairing mechanism. Additional elements are neces- 
sary and in the case of HF materials, it was suggested 
that spin fluctuations - either in the itinerant |^ or 
localized limit as magnetic excitons 0] - may replace 
the phonons as mediators of the pairing. 

In order to identify the gap symmetry and clarify its 
relationship to the pairing mechanism, selfconsistent so- 
lutions of the SC gap equations with physical momentum 
dependent pairing potentials and accurate band struc- 
tures are unavoidable (c/. [H, |^). Unfortunately, in the 
case of HF's the situation is rather complex. In addition 
to the three-dimensionality, we must deal with several 
anisotropic bands that cross the Fermi level (Ep) pro- 
ducing numerous highly anisotropic Fermi surface (FS) 
sheets. As a consequence, no realistic selfconsistent cal- 
culations of SC have been achieved so far in HF materials. 

Here we introduce a computational approach to multi- 
band SC in a HF material, by directly combining rela- 
tivistic band-structure calculations with a physically mo- 
tivated procedure to solve efficiently and accurately the 
multiband BCS gap equation selfconsistently. The effi- 
ciency of our procedure allows a systematic study of the 
influence of the momentum structure of the pairing on 
the resulting gap symmetry in an anisotropic multiband 



system. We focus on UPd2Al3 which is a fascinating HF 
superconductor having a moderately large specific heat 
coefficient 7 = 140mJ/molK^, and a relatively high crit- 
ical temperature Tc — 2K . It orders antiferromag- 
netically below Tn = 14.3 K with an ordered moment of 
O.SS/is/U-atom which is large compared to the mo- 
ments of other HF superconductors as, e.g., UPta, which 
has a moment of only 0.03 /i^ A further anoma- 

lous feature is that the antiferromagnetic (AFM) order 
coexists with SC below 2K. Coexistence of magnetism 
and SC was observed for other materials, e.g., containing 
4/ elements, but in those cases the magnetism is due to 
the localized 4/ electrons far from Ep, whereas the SC 
is carried by itinerant electrons at Ep. Conversely, for 
UPd2Al3 most of the recent studies reveal that the SC, 
magnetic order, as well as HF behavior all involve the 
uranium 5/ states Q . 

Previous investigations revealed that UPd2Al3 can be 
classified as a spin singlet, BSC-type superconductor 
p^ , pTf . The appropriate starting point is therefore the 
coupled multiband BCS gap equation p8|, for T ^ 0: 
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where A„ is the gap of band n, and i?„ is the corre- 
sponding singleparticle energy relative to Ep. Vnn' is 
the pairing potential. The coupled multiband gap equa- 
tion is at present too complicated to be solved numer- 
ically on the real FS of the material and with physical 
momentum-dependent pairing potentials. However it is a 
common experimental practice not to consider a gap sym- 
metry defined separately for each band, but to work with 
a "global" gap symmetry valid for all bands, i.e., defined 
in the whole Brillouin zone (BZ) (see, e.g., |P, 111 p3). 



Our goal is to compute selfconsistently such global gap 
symmetry. To achieve this goal, one could proceed as fol- 
lows: solve for each individual band the singleband BCS 
gap equation, and subsequently, combine and interpolate 
these singleband gaps to a generahzed gap A(fc) defined 
in the whole BZ. This is precisely what we do by solving 
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as an ansatz a modified gap equation 



A(fc) = -^ 



V{k,k')A{k') 



k' 



2[G(fc') + |A(fc') 



,1/2 



(2) 



withG(fe) ^ + + - On the FS sheet 

of band n, Gk is equal to and thus A(fc) — A„(A;). 
Note that intraband as well as interband processes are 
taken into consideration. A simplification specific to 
UPd2Al3 has been implemented in Eq. (2) which could, 
however, be relaxed. Since all bands at Ep have mainly 
5/ character , there is no need to discriminate differ- 
ent band dependencies in the pairing potential, which we 
thus assume to be identical for all n. 

The energy dispersions En{k) were computed within 
the framework of density functional theory in the local 
spin-density approximation (LSD A). Knopfle et al. 
demonstrated a very close agreement between the mea- 
sured j2^ and calculated in this framework de Haas-van 
Alphen frequencies. Consequently, the FS of UPd2Al3 
is accurately described by the LSDA energy dispersions. 
Since the relevant physics for SC happens in the vicin- 
ity of the FS, the LSDA energies provide the appropriate 
input for our investigation of the gap symmetry. We 
computed the LSDA energy bands of UPd2Al3 using the 
relativistic augmented-spherical-wave method |2^. The 
resulting band structure of UPd2Al3 along characteristic 
symmetry lines is shown in Fig. 1 for the paramagnetic 
and AFM state. The FS resulting from our AFM bands 
consists of four types of sheets, and is practically identical 
to that computed previously ||l3|, as is also the ordered 
AFM moment of 0.81 ^b, which agrees with the experi- 
mental moment of 0.85 /is ^t| . 

For the pairing potential we consider here two cate- 
gories of pairing kernels that are generally accepted and 
were shown to produce unconventional order parame- 
ters. The small-q pairing potential, adopts at small 
wavevectors k ~ k' the following form in momentum 
space V{k,k') = ■^2+(F^F)2+ - k'). This ker- 

nel is characterized by a smooth momentum cut-off 
which selects the small-g processes in the attractive part 
p3| , [2^ , p5| , p6| At larger wavevectors the repulsive 
Coulomb pseudopotential fj,*{k — k ) may prevail. This 
type of attractive interaction may be due to phonons if 
the system is close to instabilities and screening with 
short range Hubbard-like Coulomb terms is involved [ p5[ . 
It has been claimed that this type of kernel may account 
for unconventional SC in HTSC |||, |4[ ||, |26[|3 and in 
other materials, including HF's psf . Agterberg et al. |2|] 
have recently suggested for HF superconductors a mul- 
tipocket FS with an attractive intrapocket and repulsive 
interpocket potential. Such situation results naturally 
from a potential dominated by small-g attractive pairing 
as the one considered here, if the FS contains multipock- 
ets separated in momentum space. We also considered 
the case of pairing mediated by spin fluctuations using a 



phenomenological Millis-Monien-Pines pairing potential 
pof which has the form V{k, k') = ^'i_^(^j^Xk'-Q)'i where 
Q = ±(0, 0, 7r/2c) in UPd2Al3 (c is the c axis lattice con- 
stant). The Q value was verified by neutron scattering 
experiments p0| . Pairing through spin fiuctuations has 
been suggested for virtually all unconventional SC's, in- 
cluding UPd2Al3 0,0. 

To obtain a selfconsistent solution of our modified gap 
equation with such momentum dependent kernels and 
the original LSDA bands we have used a Fast-Fourier- 
Transform technique. The problem is solved iteratively 
on a symmetric part of the BZ that we discretize with a 
128 X 128 X 128 momentum grid. For each point of our 
3-D grid we have our LSDA energy bands as an input. 
Within our procedure the momentum space problem is 
fully resolved numerically without any simplification or 
bias on the resulting gap symmetry. 

The accessible even parity gap states can posses s- 
or d-wave symmetry. As UPd2Al3 and thus our LSDA 
bands obey the hexagonal Dgh point group symmetry, 
the accessible gap states of d-wave symmetry that have 
nodes should transform according to the irreducible rep- 
resentations of the Dgh group |^ shown schematically 
in Fig. 2. In our calculations we used the small-g phonon 
pairing potential with various parameter choices for both 
the paramagnetic and AFM energy bands and the spin- 
fluctuations pairing potential only for the AFM bands. 
An initial gap configuration is chosen randomly and the 
system converges within the iteration cycle towards the 
most favorable solution. Also, we adopted as initial gap 
each of the representations in Fig. 2 and tried if the sys- 
tem would converge towards the chosen configuration. 
Depending on the parameters in pairing kernel {q^, V, 
fi*), we can achieve as selfconsistent solution either a gap 
having s-wave or d-wave symmetry. Quite unexpectedly, 
in the latter case we find that the only accessible gap state 
with nodes in UPd2Al^ belongs to the irreducible repre- 
sentation Aig. No other representation possessing nodes 
was accessible for all the parameters in both types of 
pairing kernels. 

We show in Fig. 3 some examples of our selfconsis- 
tently computed gap functions in the plane H — A — H, 
H — K — H which contains the z axis (the A — F — A 
axis) and in the plane obtained by a 7r/6 or 7r/2 rotation 
around the z axis (the plane L — A — L, L — M — L). All 
selfconsistent solutions shown in Fig. 3 have two lines of 
nodes perpendicular to the A — F — A axis, thus belonging 
in the Aig representation. To study the influence of the 
coexisting AFM order on SC we have made calculations 
using both paramagnetic and AFM bands. In both para- 
magnetic and AFM cases, the only solutions with gap 
nodes are of Aig type. We note that our results with 
the spin-fiuctuations kernel confirm the model analyses 
performed by Huth et al. |l^] and Bernhoeft [Q. In 
particular, Bernhoeft predicted for SC mediated by spin 
fluctuations in UPd2Al3 a characteristic symmetry prop- 
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erty of the global gap, viz. A(fc) = —A{k + Q), which, in 
its simplest form, would be fulfilled by a gap having the 
functional dependence A(A;) oc cos(A;2c). Our gap com- 
puted selfconsistently with the spin-fluctuation pairing 
kernel does obey this symmetry property, and the shape 
of the gap does correspond closely to cos(fczc) (see Fig. 3). 
In addition, we find here that small-g phonon pairing 
may lead to the same d-wave Aig state, demonstrating 
that there is not a one-to-one correspondence between the 
gap symmetry and a specific microscopic mechanism, al- 
though experiments suggest that spin degrees of freedom 
are likely involved in the pairing ^, |3^, |3^ . 

There is overwhelming experimental evidence that in- 
deed UPd2Al3 has a gap with the Aig node structure. 
Recent tunnel measurements along the z axis showed the 
absence of a node in this direction However, the 

presence of nodes is definitely established for UPd2Al3, 
and since Aig is the only even parity gap with nodes 
which is nodeless along the z direction (see Fig. 2), it 
was deduced a posteriori that the gap symmetry must 
be Aig type In addition, the observation |^, |3|] 

of a spin-fluctuations peak below Tj, at (0, 0, 7r/2c) would 
be forbidden by the involved coherence factors unless the 
gap changes sign along the z axis as in the Aig represen- 
tation. Finally, measurements of the angular dependence 
of the critical field indicate two lines of nodes perpendic- 
ular to the z axis |]l9| again suggesting the Aig represen- 
tation. UPd2Al3 is one of the few HF superconductor for 
which such a clear picture of the nodal structure of the 
gap is obtained from the experiments. We consider the 
agreement of our results with the experiments as strong 
support of our calculations. 

The surprising robustness of the Aig solution results 
from the following general rule: If a system must choose 
a gap symmetry with nodes because of the repulsive ef- 
fective interaction at large wavevectors (short distances) 
it chooses the representation that has the minimum num- 
ber of nodes as far as possible from the high-density ar- 
eas. In fact, the system must maximize the condensa- 
tion free energy and this is obtained when there is a gap 
in the high-density areas and the node (gapless) areas 
are kept minimal. Examining the paramagnetic bands 
of UPd2Al3 (Fig. lb) we see that the high-density areas 
near the FS due to saddle points in the band dispersions 
are found essentially near the A point and near the H 
point. The Aig representation prevails because it is the 
only one without a node near the A point. In the case of 
the AFM bands (Fig. Ic), the saddle points of the bands 
near the FS are found essentially in the z = plane (near 
the r point, along the F — K, K — M and M — F symme- 
try lines) and near the F point along F — A' (A' = A/2). 
The high-density areas in the z = plane exclude the 
Big, Eig and B2g representations which have a line node 
in the z — Q plane. From the remaining Aig, E2g and 
A2g representations, Aig is the only one without a line of 
nodes parallel to the A — F — A axis which would cross the 



high-density areas at F and near F along F — A'. Also, 
due to the AFM symmetry a high density occurs again 
at A, as A is equivalent to F. The node perpendicular to 
the F — A path in the Aig representation is close to A' 
and therefore does not cross the high-density areas near F 
and A. Note that A2g, Big and i?2g are also handicapped 
by the fact that they have more nodal areas. 

To stress the generality of this argument we recon- 
sidered the case of cuprate HTSC. Saddle points pro- 
duce high-density areas in the vicinity of (0,7r/a) and 
symmetry related areas (a is the basal-plane lattice con- 
stant). The dr^2_y2 gap symmetry corresponds to nodes 
along the (tt/o, tt/o) direction as far as possible from the 
high-density (0, ±7r/a) and (±7r/a,0) areas. Using both 
types of pairing kernels as adopted for UPd2Al3 (now 
Q = (tt/o, 7r/a)) we were never able to obtain a dxy solu- 
tion for which the nodes cross the high-density (0, tt/o) 
areas no matter the details of the interaction. 

In conclusion, we have proposed an efficient compu- 
tational approach to the multiband gap equation. Our 
selfconsistent calculations of the gap in UPd2Al3 demon- 
strate that the only accessible gap symmetry with nodes 
transforms according to the Aig irreducible representa- 
tion of the Dg/i point group, independent of whether the 
pairing potential is dominated by small-g attractive pro- 
cesses involving phonons or by spin fluctuations. The 
robustness of the Aig representation is analyzed to be 
due to the presence of high-density regions of the rele- 
vant bands in the vicinity of the F and A points, that 
preclude node formation at these points. As a general 
rule we obtain that nodes must be as far as possible from 
high-density areas in the phase space. 
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M. Huth, M. Lang, P. Thalmeier, N. Sato and F. Steglich 
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figure of the irreducible representations of the Dgh group. 
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FIG. 1: (a) The hexagonal Brillouin zone with high symme- 
try points. In the AFM phase the BZ is reduced by a factor 
two along the z axis (denoted by the primed letters), (b) The 
energy bands along high symmetry directions in the param- 
agnetic, and (c) the AFM phase of UPd2Al3. The relevant 
bands for superconductivity are shown by the dotted lines. 



FIG. 2: The even parity d-wave-type irreducible represen- 
tations of the D(,h point group. While in principle all are 
accessible for the gap, only the A\g results from our selfcon- 
sistent gap calculations. 



FIG. 3: Examples of the selfconsistent solutions of the gap 
equation in characteristic planes (with sides indicated) con- 
taining the r point (in the center) and the z axis. The thick 
black lines show the nodes and in the dark gray areas the 
gap is negative. All nodes cut the z axis as in the A\g repre- 
sentation. Shown are gaps computed with the small-g pairing 
and paramagnetic energy bands (a), and the spin-fluctuations 
pairing and AFM bands, (b) and (c). 
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